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Abstract 

We investigate the exact results for circular 1/4 and 1/2 BPS Wilson loops in the d = 3 
Af = 4 super Chern-Simons-matter theory that could be obtained by orbifolding Aharony-Bergman- 
Jafferis-Maldacena (ABJM) theory. The partition function of the Af = 4 orbifold ABJM theory has 
been computed previously in the literature. In this paper, we re-derive it using a slightly different 
method. We calculate the vacuum expectation values of the circular 1/4 BPS Wilson loops in 
fundamental representation and of circular 1/2 BPS Wilson loops in arbitrary representations. We 
use both the saddle point approach and Fermi gas approach. The results for Wilson loops are in 
accord with the available gravity results. 
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1 Introduction 


In the AdSs/CFT 4 correspondence [l]~3], 1/2 BPS Wilson loops in d = 4 J\f = 4 SU(N) super Yang- 
Mills theory are dual to fundamental strings in type IIB string theory in AdSsxS' 5 spacetime [4 ]-[t] . 
When string theory is weakly coupled and the supergravity approximation is a good one, the dual 
d = 4 AT = 4 super Yang-Mills theory is strongly coupled. To compare with the gravity results, one 
has to know the vacuum expectation values of the Wilson loops at strong coupling. To do this it was 
proposed in [8,9| that d = 4 AT = 4 super Yang-Mills theory is related to the Gaussian matrix model, 
and this was proved in [lOj using localization techniques. 

There is a similar but more complicated story in the AdS 4 /CFT 3 correspondence. M-theory in 
AdS 4 xS 7 /Zfc spacetime, or type IIA string theory in AdS 4 xCP 3 spacetime, is dual to the d = 3 N = 6 
super Chern-Simons-matter (SCSM) theory with gauge group U(N ) x U(N) and levels (A:, —k), which 
is known as Aharony-Bergman-Jafferis-Maldacena (ABJM) theory 111. In ABJM theory there are 


1/6 BPS 12-14 and 1/2 BPS [151 Wilson loops. The 1/6 BPS Wilson loops are closely related to the 
1/2 BPS Wilson loops in AT = 2 SCSM theory in 116]. Localization techniques have been applied to 


ABJM theory and other SCSM theories with fewer supersynmreties 17-19 and lead to matrix models 
that are more complicated than the Gaussian matrix model. By using localization, one can calculate 
the partition function and vacuum expectation values of Wilson loops at both weak coupling and 


strong coupling 15,17,20-22 . The computations in 22 are based on the saddle point solution of the 
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ABJM matrix model at large N limit with finite k, and we will call such a method the saddle point 
approach. Furthermore, the ABJM matrix model could be reformulated as an ideal Fermi gas with 


a complicated potential 23 , and one can calculate the vacuum expectation values of BPS Wilson 


loops with fixed winding number using the Fermi gas approach |24| . One can also use the Fermi 
gas approach to calculate the vacuum expectation values of the 1/2 BPS Wilson loops in arbitrary 
representations [25]. 

By Z r orbifolding the U(rN) x U(rN) ABJM theory, one can obtain a d = 3 Af = 4 SCSM 
theory with gauge group JJ(N) 2r and levels (k,-k,--- ,k,—k). This theory is dual to M-theory in 
AdS 4 xS 7 /(Z r xZ r fc) spacetime 26-29 . The partition function of the orbifold ABJM theory has been 


calculated using Fermi gas approach in 1301, and in this paper we will re-derive it using a slightly 
different way. In the orbifold ABJM theory there are 1/4 and 1/2 BPS Wilson loops, and the 1/2 
BPS Wilson loops in fundamental representation should be dual to M2-branes with one dimension 


wrapping on the M-theory circle 31,32 . In this paper, we will calculate the leading contributions of 
vacuum expectation values of the Wilson loops using the saddle point approach in the large N limit 
with k and r being finite. We will also calculate the perturbative pari]/] of the vacuum expectation 
values of the Wilson loops using the Fermi gas approach. The results are in agreement with the 
available gravity results. 

In the Af = 4 orbifold ABJM theory with gauge group U(N) 2r , there are 2 r linearly independent 
1/2 BPS Wilson loops that preserve the same supersymmetries, but there are not so many 1/2 BPS 
branes in M-theory in AdS 4 xS 7 /(Z r xZ r fc) spacetime. It was conjectured that these Wilson loops are 
1/2 BPS classically, and only a special linear combination of them is 1/2 BPS quantum mechanically 
j|321. If all the 2r Wilson loops are 1/2 BPS and each of them differs from 1/4 BPS Wilson loop by 
a Q-exact term quantum mechanically, we can calculate their vacuum expectation values in a matrix 
model as shown in this paper. If only a special linear combination of the 2 r Wilson loops is 1 /2 BPS 
quantum mechanically and it differs from an 1/4 BPS Wilson loop by a Q-exact term, we can calculate 
its vacuum expectation value in the matrix model. If there is no 1/4 BPS Wilson loop that differs 
from such an 1/2 BPS Wilson loop by a Q-exact term, we cannot calculate the 1/2 BPS Wilson loop’s 
vacuum expectation value using currently available localization techniques. In this case, the large part 
of the calculation in the paper are just some results in the matrix model and have nothing to do with 
vacuum expectations values of half-BPS Wilson loops. 

The rest of the paper is arranged as follows. In Section 2 we review the results in ABJM theory, 
including the partition function and vacuum expectation values of Wilson loops. In Section 3 we 
investigate the partition function of the AT = 4 orbifold ABJM theory. In Section 4 we review the 
circular 1/4 and 1/2 BPS Wilson loops in the Af = 4 orbifold ABJM theory in Euclidean space. In 
Section 5 we calculate vacuum expectation values of Wilson loops with fixed winding number using 
the saddle point approach. In Section 6 we calculate vacuum expectation values of Wilson loops in 
arbitrary representations using the Fermi gas approach. We end with conclusions and discussions in 


: By this we mean to include all of the 1/N corrections, putting aside the non-perturbative contributions. 
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Section 7. In Appendix A we investigate if there are more general 1/2 BPS Wilson loops in A/" = 4 
orbifold ABJM theory other than the ones found in 131,321. We find no new ones. 


2 Results in ABJM theory 

In this section we review some results in ABJM theory. This includes the partition function and 
vacuum expectation values of circular 1/6 and 1/2 BPS Wilson loops. We focus on what will be used 
in the following sections, so this is merely a brief review. 


2.1 Partition function 

The partition function of ABJM theory with gauge group U(N) x U(N ) and levels (k,—k) can be 


localized to be the ABJM matrix model 17 


N 


= w i n 


d fii d Vi 11 i<j 


; ru 2 sinh 




2 sinh 


2=1 


2n 2n 


n ,(2cosh^ 


x exp 


ik 

47T 


I] (Pi - ^) 


( 2 . 1 ) 


The partition function of the matrix model (2.1) can be written as the canonical partition function 


Z(N) of 7V-particle free Fermi gas with the one-particle density matrix being 23 

P = e ~^ j 


( 2 . 2 ) 


whose explicit form will not be used in this paper. Note that H is the one-particle Hamiltonian 
operator. To calculate Z(N), one can firstly calculate the grand partition function 

+oo 

S(aO = X] z n Z(N), (2.3) 

N=0 

with Z(0) = 1, z = e^ being the fugacity and fi being the chemical potential. The grand potential is 
defined as 

J(p) = logS(ju). 

And then one gets 


Z(N) = r 
V ’ 2vri 


(2.4) 

(2.5) 


One can define j{p) according to [33 


e J M = 


E' 

/= — OO 




and then 






One adopts the phase space formulation of quantum mechanics, and defines 

n (p)= j ^e(p-H)w, 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 
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with h = 2nk, 9{x) being the Heaviside step function, and W being the Wigner transformation. 
The quantity n(fi) counts the number of one-particle states whose energy is less than /i. Using 
the Sommerfeld expansion one can get the expectation value of particle number N(/j) in the grand 
canonical ensemble 

N(/i) = Trd/j, csc('Kd fi )n(fi). (2.9) 

It is standard in the grand canonical ensemble that 

dJ{n) 


and so we get 


We find that when /i -» — oo, 


"<"»= a, • 

J(/i) = f N(u)du. 

J — OO 

0, i = 0,1,2,- 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


This results is very useful for us. Note that the way from n(/r) to N(/j) and then to J(//) is a slightly 
different method of getting J(/r) to the one in the original paper |23j. 

In the large // (i.e. large N) limit, one can split a quantity into the perturbative part and non- 
perturbative part. The perturbative part is denoted as pt. The non-perturbative part is exponentially 
suppressed in the large /j (i.e. large N) limit, and it is denoted as np. In this paper we will mainly 
focus on the perturbative part. It turns out that [23 


n pt (/x) = C/i 2 + n 0 , 


with 


One then gets 


C = 


-K 2 k ’ 


n °“~3fc + 24- 


JVpt(/i) = C[x 2 + B, 

C , 

Jpt(v) = +B[i + A, 


(2.13) 


(2.14) 


(2.15) 


where 


B = n 0 + 


tt 2 C 1 k 
“ 3k + 24' 


(2.16) 


Here A appears as an integral constant, and its exact form depends on the full form of N(fj,). One 


can find the result for A in 134,35 . One has 


jpt(k-) — 


and then one gets the perturbative part of the partition function 123,36 

Z p t(JV) = U- 1/3 e^Ai[C' 1/3 (lV - B)], 
with Ai(x) being the Airy function. 


(2.17) 


(2.18) 
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2.2 Wilson loops 


The representations of group U(N) and supergroup U(N\N) can be denoted by Young diagrams. We 
write a general Young diagram as R, and it can be a representation of U(N ) or U(N\N). 

We consider the hook representation R = (a|i>) with a + 1 boxes in the first row and one box in 
each of the remaining b rows. For both the 1/6 BPS and 1/2 BPS cases, a Wilson loop with winding 
number n is related to Wilson loops in the hook representations by 


71—1 


W n = ^2(-i) b w { 


(n—l—b\b) ■ 


(2.19) 


6=0 


When n = 1, it is just the fundamental representation. In the matrix model (2.1), the circular 1/6 
and 1/2 BPS Wilson loops with winding number n can be written as |15,17 


<Wr /6 ) = ( £ e""), <W? /6 ) = ( £ e nVi ), (W? /2 ) = ( £ [e nw - (~) n e w *] ), (2.20) 


with n being the winding number of the loop and the right hand sides being the expectation values 
in the matrix model. For their expectation values one has the relation 


<wT/ 2 > = <"?/6> - (-mr /6 > = (wr /6 ) - (-mr /6 >*, 


( 2 . 21 ) 


with * being the complex conjugate. 

In the large N limit with hnite k, i.e. the M-theory limit, the values pi and ui at the saddle point 


can be denoted as a continuous distribution 22 


/j,(x) = VNx + u{x) = VNx — 

4-7T 47r 


( 2 . 22 ) 


with the uniform density 


p{x) = 


2x* 


x G [—a;*,a;*], x* = 7r 


(2.23) 


In the saddle point approach the Wilson loop vacuum expectation values can be calculated as 




1/6 / 
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0 n[i(x 


^ p{x) d, 
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' — X * 

rx* 


2 mr V 2 


X nnV2\ 


—e 




i/2> « N 


i n 1 k 


& n^x) _ ,_™ e nv{x) ^ i_// e n^Y2A_ 

J 4n7r 


(2.24) 


The exponentially suppressed terms are omitted here. Note that one can only get the correct leading 
contribution of large N in the saddle point approach. 

The vacuum expectation values of circular Wilson loops can also be calculated in the Fermi gas 


approach 24 . One firstly calculates 

m(jn) = 


d< 7 d 7V n( V p) 

—0(n~H) w e * , 


(2.25) 


and then using Sommerfeld expansion one gets the 1/6 BPS Wilson loop expectation value in the 
grand canonical ensemble 

M(p) = 7t<9 m csc(7 Td fl )m(p). (2.26) 
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Then the 1/6 BPS Wilson loop expectation value with winding n in the canonical ensemble is 




Similar to the partition function, one has 




1 


dh e J pt (M)-^ Mpt ( / ,) ; 


1/6/pt Zpt(N) i„ ioo 27ri" 
with non-perturbative contributions being neglected. It turns out that 


2 n/i 


m p t (p) = (Dp + E)e 

2irn 2irn 2 n M 
M pt (p) = esc —r-e k 


with 
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k 2nn . „ _ 

p + — - vr cot —) A + A 
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i n+1 k 
47r 2 n 
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- - i// r 
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Here H n is the harmonic number, 

n 

Hn = ^d : 

d= 1 

with I?o being 1. Then the 1/6 BPS Wilson loop vacuum expectation value is 

A 2 l- 1 / 3 (M _ A _ 6n+l \ 

ViPfc/ P v 24 3fc ) 

. _ 1 _ 

7T 2 fc 


/ o \ —i/3 Ai' 
<"?/6>Pt = - I — I F 


Ai 




Ai 


+G- 


( 2 j~ 1//3 lV 6n+l \ 

ViPfc/ P v 24 3fc ) 


Ai 




where 


27TO 27TO 

A = —-— esc- D, 

k k 

2irn 2 ^rn 

G = —— esc 


k k 

The 1/2 BPS Wilson loop vacuum expectation value is 


k 2nn , „ „ 

)D + E 


(W, 


1/2)pt 


•n—1 


2 ttti 


Ai 


CSC 


f/V _ A _ 6n+l 1 
\7v^k) V JV 24 3fc / 


k 


Ai 


(AWV-A-s) 


(2.27) 


(2.28) 


(2.29) 

(2.30) 

(2.31) 


(2.32) 


(2.33) 


(2.34) 


Now we turn to Wilson loops in hook representations based on [25]. There the density matrix for 
Fermi gas dual to ABJM theory was obtained as 

1 „ 1 


p = s/QP\/Q , with P = 


2 cosh | 


, Q = 


2 cosh | 


(2.35) 


Though it is the same as the matrix in 123 but different from the one in 24 , it gives the same partition 


functions and vacuum expectation values of BPS Wilson loops. One of the key steps in 25 is the 
following result 

/( e wivGC 


sM(n$P> + 


(2.36) 
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where 


pf = VQ f^ w ^ p f( w )VQi with w = e ^ P 


The density matrix pf with f(W) = (1 + tW)/(l — sW ) can be written as 25 

OO 

p f = p+(s + t) Y s a t b \b)(a \, 


(2.37) 


(2.38) 


a, 6=0 


where |a) and (6| are defined in the coordinate q representation as 


(qW) = 




2 cosh | 


, (b\q) = (q\b)* = 


e ( b +|)f+T b (*>+i) 


yj 2 cosh | 


(2.39) 


For the half BPS Wilson loop in a hook representation (a|7), the generating function is given 
by |37| 


°o JJ N , 

i+(*+*) e = sdet (rwrff)= n ?! 


(1 + te^)( 1 + se^) 


a.b =0 


^ (1 — se^)(l — te v i)' 


(2.40) 


with 7/ = diag(C7 /i , — U u ), U M = diag(e Mi ), U v = diag(e^). Therefore, the grand canonical ensemble 
expectation value of 1/2 BPS Wilson loop generating function in ABJM theory becomes 


1 + (s + t) sCl t b W(a\b) 

a,b =0 


GC det(l + zpf) 


det ( 1 + (s + t) Y^ sa t b 


a,b =0 


det(l + zp) 

OO 

2 |(,){a|)=l + ( S + t)^ 8 “t l >{o| I 


1 + zp 


a,6=0 


+ zp 


I b). (2.41) 


One gets the relation 




1+ zp 


|5><a| =Tr 


e H \b)( a \). 


+ 1 


(2.42) 


As discussed in |25], the perturbative part of the half BPS hook Wilson loop in ABJM theory is 
determined by the topological vertex of C 3 in 


38 


tw \GC = q^ a{a+l) ^ b(b+1) a+b+l 2( a+ 6+iV 
( W {a\b)) pt [ a + ft+1 ][ 0 ]![ 6 ]l 


(2.43) 


TX r n 1 L _ ]± 

with q = e k and [n\ = q 2 — q 2 . 

Let us consider the circular half BPS Wilson loops in non-hook representations. One can decompose 
the Young diagram for a non-hook representation into hooks from the upper left to the lower right to 
get (ai|&i), • • • , (o s |6 s ). This general representation will be denoted as R = (01 • • • o s |6i • • • b s ). The 
Giambelli formula states that 


bL( aia2 ... a;s |fe 1 b 2 ...b s )(e^% e J ) detW( ap , e j ). 


(2.44) 


2 In the remainder of this section, we will only discuss the half BPS Wilson loop and omit 1/2 in the subscript. 




















The authors of 25 considered the following generating function 


The computations in [25] give 


W(N) = (det[5 pq + tW {aplbq) (e^,e^)}). 


W{N) = m / ^\[dvi]Detp f {y u Vj), 


(2.45) 


(2.46) 


where 


d Hi fikp. 


[d/x*] = ex P 


47T 


’ = ^ exp ( ~~ 


iku[ 
47T J ' 


Pfiyuvj) = / [M 


_|_ + \ p(b P +l/2)i/i (a p +l/2)fi \ 

2 c 0 sh^ V 2 cosh^ ' J' 


p =i 


Then one has 


z n W(N) = Det(l + zp f ). 


(2.47) 


(2.48) 


N=0 


The multiplication between boldface variables is understood as matrix multiplication with indices p, v 
and summation being replaced by integration with measures [d/x], [dz/]. Then by introducing 


((a|-T)w=a ( - 1/ ^.(-^| ( ,))M = ^ 1/ ^. 


(2.49) 


one can get 


Det 


s 1 

1 + z(^P + -^q\ b p)( 

— 1 


\TQ 


Q 


= Det(l + zp) det \5 vq + zf(a„|(l + zp) 1 16„)1. (2.50) 

pq l rH r h j 


Taking t = 0 in the above results, one gets 


(z) = Det(l + zp), 


(2.51) 


and then one has 

(det(«Jpg + tW {ap \ bq) )) GC = det[(5 P q + zt(a p \(l + zp)~ 1 \b q )}. (2.52) 

The coefficient of t s in both sides of the above equation gives 

{W {aia2 ... as \ blb2 ... bs) ) GQ = det [z(a p \(l + zp)~ 1 \b q )]. (2.53) 

PiQ 

Restricted to hook representation cases, one has 

(W {alb] ) GC = z(a\(l + zp)- 1 \b). (2.54) 

So finally one gets 

(^(a 1 a 2 -a s |6i6 2 -b s )) GC = de t( W (a P |6 9 )) GC - ( 2 - 55 ) 

This shows that, for the non-hook representation cases, the expectation values of half BPS Wilson 
loops in the grand canonical ensemble can be written as the determinant of the expectation values of 
Wilson loops in the hook representations. In other words, they are Giambelli compatible. 
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3 Partition function 


The computation of the partition function of the J\f = 4 orbifold ABJM theory can be localized to the 
matrix model |17j^ 


. r—1 N 


Z r {N) = 


N\ 2r 


n n dm, n IL<,• (2 )' ( 2 ' 

1111 2vr 2vr H 


t =o i=l 


T 11, ; (2cosh ( 2 cosh Vt+1 ' i 2 


x exp 


ifc 

47T 


I] (pli - 


i.i 


(3.1) 


When r = 1 it is reduced to the ABJM matrix model (2.1). It can be written as canonical ensemble 

(3.2) 


partition function of an ./V-particle Fermi gas with one-particle density matrix 123 

-rH 


Pr — 6 


with H being the same as that of ABJM theory in (2.2). 

We calculate the partition function in the Fermi gas approach. We firstly have 

dqdp 


n r {p) = 


27 xh 


9{p - rH) w 


(3.3) 


dqdp ( p 


27 xh 


e(^-H) = n (—) 

V r / w V r / 


with n(p) being the same function as (2.8). Then 


N r (p) = 7r9 M csc(7r d^)n r {p) = 


sin(rTr^) N /p 
r sin(7r5 A1 ) \rJ 


Then using (2.12) we can get 


J r (p) = / N r (u)du = 

J — OO 


sin(r7r9 M ) /p 


sin(7r9u) V r 


Note that we have the following expansion 


sin(r 7 r 3 At ) - 27^.2 


r sin^c^ 
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+ 


360 


(3.4) 


(3.5) 


(3.6) 


Formula (3.5) is a convenient way to get the grand potential J r (p) with Hamiltonian rH from the 
grand potential J(p) with Hamiltonian H, including both the perturbative and non-perturbative parts. 
From the results of ABJM theory we have 


n, 


?\p) = c r p 2 + n r 0 , N?\p) = C r p 2 + B r , 


C r 


B r p + A r , 


with 


«. = « = 


r r 2 TT 2 r 2 k 

ir 2 C(r 2 - 1) 


’ n ° = n ° = -3fe + 24- 


(3.7) 


(3.8) 


B r = B — 


r 2 — 2 k 


3 r 2 


3 r 2 k 24 


J- ~ ; A r — r A . 


3 Supersymmetric localization in d = 3 JV = 2 SCSM theories was first studied in 
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Then we have the perturbative part of the partition function 


Zf(JV) = (£) 


- 1/3 


rA 


Ai 


- 1/3 


(3.9) 


This is in accordance with the result in |30J, and here we re-derive it in a different way. 

The non-perturbative part of the grand potential for ABJM theory J np (/i) is a summation of terms 
of the form [23] 

(a // 2 + 6/1 + c)e _d/i , (3.10) 

with a, b, c, d being constants and d > 0. Correspondingly in the grand potential of the AT = 4 SCSM 
theory J/ ip (/i) there is the term 


sin(r7n9 M ) / a 2 & 

sin(7r<9 A 


y + e ^ -A + br/i + c r ) e 


(3.11) 


with 


a r — 9 /r I 


Or - z, i -) - 5/: a 


b , f d \ a „ 


Cr = cfr(z)-zfr(z)+ ~2 fr ( Z 


Here we have defined the function 


, , \ sm(r7rx) 
Jr[x) = — 


sin(7rx) 


(3.12) 

(3.13) 


Note that when x = l is an integer we have 

fr{l) = r(-) (r - 1)Z , m = 0, 

m = ~(-) ir - 1)l (r 2 - 1) 


(3.14) 


4 Circular BPS Wilson loops 

In this section we review the circular 1/4 and 1/2 BPS Wilson loops for the AT = 4 orbifold ABJM 


theory in Euclidean space [31,32 . This theory is an SCSM theory with gauge groups U(N) /r and 
levels (k, —k, ■ ■ ■ , k, — fc)j^] In d = 3 Euclidean space we use the convention of spinors in [39], and 
especially we have the coordinates x^ = (x 1 ,® 2 ,^ 3 ) and the gamma matrices 


'fj = (--W 3 ), 


(4.1) 


with cr 1 ’ 2 ’ 3 being the Pauli matrices. The circle is parameterized as x^ = (cost, sinr,0). 

Using every gauge field A/ ; with i = 0,1, ■ ■ • ,r — 1 and matter that couples to it, one can 
define the 1/4 BPS Wilson loop 

ll// 11 ' = TrPexp (-i/drp( ( 2 ,+ 1 >(r)') , 

Z 2,+1) = + y (m- *? +1 % +1) + 


M 1 - = M l - = diag(i, -i). 


(4.2) 


4 Without loss of generality, we assume k to be positive. 
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The conserved supersymmetries can be denoted as 


v n = i l 3 e u , v 22 = -i l 3 e 22 , 
e 12 = e 2i = v 12 = P 2i = o, 


(4.3) 


where the spinors 0 u and i3 u with i = 1,2 and i = 1,2 denote the parameters of Poincare and conformal 
supersymmetries, respectively. 

Also using every gauge field A/, ; with i = 0,1, • ■ ■ , r — 1 and matter that couples to it, one can 
define the 1/4 BPS Wilson loop 


wffi = TrV exp 


-i j) dri^ 2 ^(r)j , 


2vr 
~k 

N/ = N/ = diag(i, -i). 


(4.4) 


This kind of 1/4 BPS Wilson loop preserves the same supersymmetries as the previous one (4.3). The 


1/4 BPS Wilson loops (4.2) and (4.4) can be combined to give a 1/4 BPS Wilson loop 


W i/\ = TPPexp ( -i <p d tL\J a (t 


L m _ ( 

-1 


J l/4' 


(4.5) 


Using two adjacent gauge groups A/ . A/ ; in the quiver diagram and matter that couples to 
them, one can define the 1/2 BPS Wilson loop 


W[% = AAV exp ^ j dTlif / ) 2 (r)'j , 
[(2£+l) fW) \ 

f m 1. 


At) _( 

h l/2 ~ 


A^ + 1) = Af +1 )^ + y (M^f +1) ^' 


(2£+l) + M j&i ^(2 £)) 1^1 


a (2() = a^x^ + y (^>(2^1)^ 1} + at-v; 2t)^f £) ) l®l. 
/i (2£) = y^7 ( 2 *V(2#l> = \J^r$ i 2 \ 2 *)N, 


(4.6) 


= dia §(h —i)> M j = = diag(-i, -i), 

^ )q = ^(e iT/2 , e _ir/2 ), r, (2 , )Q = (e"^ 2 , e 1 ^ 2 )/?, P/3 = i. 

Note that f3 and (3 are Grassmann even constants. The conserved supersymmetries are 

= Z730 1 *', i? 2i ~ = —i^O 21 , i = i, 2. (4.7) 

In the terminology of |32|, the above 1/2 BPS Wilson loop is called the -i/q-loop, since it is coupled to 
fields V ; ( 2 £) and ■ 
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Similarly there is a 1/2 BPS 02-loop that is coupled to fields '4^2f.) an< ^ ^2 


m 


of 31 such a 1/2 BPS Wilson loop can be constructed as 


Wfj 2 = TrV exp ( -i 


drL$ 0 (r) 


T/2' 


f oo _ 

^1/2 _ 


_ # +1) ft 


m 


f, 


(2 i) 


A (2£) 

2tt 


#*> = if ^ + f + N^\ 2e /; e) 


fV = 







( 2 £)\ 


A2£) _ 
J 2 ~ 


2t) I®I, 


M l j = N/ = diag(i, -i), 
fjW* = p(e iT / 2 ,-e~ iT / 2 ), 

M = i- 


T" 

M i = A? { J = diag(i, i), 

= (e" ir/ 2 ) -e-/2 ) /3, 


3 

r ?(2f)c 


32 


In the conventions 


^(2m) = Af+ 1 )^ + ^ (M^f' fl) 4, +1) + 


(4.8) 


The 02-loop preserves the same supersymmetries as the 0i-loop (4.7). 

It has been checked that the difference of 1/4 and 1/2 BPS Wilson loops is Q-exact with Q being 


some supercharge preserved by both the 1/4 and 1/2 BPS Wilson loops [ 31,32 . This applies to both 


the 0i-loop (4.6) and -02-loop (4.8), and explicitly one has 

w[% - w}% = QV( £ \ w\% - w\ l l = QV&, 


1/4 


1/2 


1/4 


(4.9) 

with and being some operators. It was conjectured that the 2 r Wilson loops W/f an ^ iff 
with £ = 0,1, • • • , r — 1 are 1/2 BPS classically, and only a special linear combination of them is 1/2 
BPS quantum mechanically [32). If it is the case, we may denote such a true 1/2 BPS Wilson loop as 

r— 1 

(4.10) 


»T/” = E («<( + 




£=0 


with eg and eg being some to-be-determined constants. Here superscript qm means that the Wilson 
loop is 1/2 BPS quantum mechanically. We do not know if Wilson loops (|4.2[), (4.4) and (4.5) are still 


(4.11) 


BPS quantum mechanically, but we expect that at least there is the 1/4 BPS Wilson loop 


r— 1 


wT 4 = ^2(c e + c e )W^ 


1/4- 


£=0 


In this case, equations (4.9) would also be spoiled. We expect that 

r— 1 


W\%-W$ = QV, with V = ^(ctV^ + ce 


(4.12) 


1=0 


Now we have three possibilities. The first is that Wilson loops (4.6) and (4.8) are 1/2 BPS quantum 


mechanically, and equations (4.9) also hold quantum mechanically. In this case the Wilson loops 


(4.6) and (4.8) have the same vacuum expectation values. The second possibility is that only Wilson 


loop (4.10) is 1/2 BPS quantum mechanically, and (4.12) holds quantum mechanically. The third 


possibility is that Wilson loop (4.10) is 1/2 BPS quantum mechanically, but (4.12) does not hold. 
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5 Wilson loops in saddle point approach 


In this section, we compute the vacuum expectation values of Wilson loops with fixed winding number 
based on the saddle point approach. 


If equation (4.9) holds quantum mechanically, we have the relations between vacuum expectation 


values of Wilson loops and expectation values in the matrix model (3.1) 




>ni/£. 


wl% n ) = w[% n ) = (E ^ ), 


(5.1) 


with n being the winding number. From Z r symmetry of the matrix model (3.1), we have 




<<•”> = <<!”> = (£ 


= n/i 0 ,i _ ^_^rigTO'o.ij \ 


(5.2) 


If quantum mechanically we have (4.12), we get 


(W^ n ) - (-Ye™ 0 ’'} y 


with 


r —1 


C = 


E ( c ^+• 


(5.3) 


(5.4) 


£=0 


It is possible that (4.12) is not true quantum mechanically. But it is still an interesting problem in its 


own right to calculate the expectation values in the matrix model 

</ Y e n »°A, / J2 A* 10 '' - (-) n e ni/ °’ 1 ] )• (5.5) 

i i 

We calculate the Wilson loop expectation values in the saddle point approach. For the matrix 


model (3.1), at the saddle point we have 22 


M®) = M®): U Y X ) = v ( x ). 


(5.6) 


with n(x) and v(x) being the same as (2.22). Then for the 1/4 and 1/2 BPS Wilson loops there are 
leading contributions of the vacuum expectation values which are the same as in the ABJM case 




(5.7) 


W) n \ _ /TT/W’ n \ ^ rnrV2\ 


(K/f) = WW) 


1/2 


4n7r 


Note that there is no i or r dependence in this result. If the Wilson loops (4.6) and (4.8) are not 1/2 
BPS quantum mechanically, we cannot use the matrix model to calculate their vacuum expectation 
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values 32 . If (4.12) holds quantum mechanically we can repeat the above process for the true 1/2 


BPS Wilson loops (4.10), and we have 




\n —1 


ck 


t\/2A 


4n7r 


(5.8) 


with constant c being (5.4). The Wilson loops in the fundamental representation are those with 


winding number n = 1. In 31 32 it was shown that a suitably positioned M2-brane in M-theory in 
AdS 4 xS 7 /(Z r x Z r fc) spacetime can be 1/2 BPS, and for the regularized on-shell action of the M2-brane 
in Euclidean space one has 


,—Sm 2 


t\/2A 


(5.9) 


We find matches of the matrix model and gravity results. 


If both (4.9) and (4.12) are spoiled by quantum corrections, the matrix model calculations here 
would have nothing to do with vacuum expectations values of Wilson loops. 


6 Wilson loops in Fermi gas approach 

In this section we use the Fermi gas approach, and study vacuum expectation values of Wilson loops 
with fixed winding number and of 1/2 BPS Wilson loops in both hook and non-hook representations. 


6.1 Wilson loops with fixed winding number 

We calculate the Wilson loops expectation values in Fermi gas approach. We firstly calculate 

mrifi) = J 

And then we can get 


d<?dp n(g + p) [ d( l d P Q (P tt\ "(g+P>. (P 

‘ = h^h > (r- H )w e =m \r 


171; 


P \p) = ( -n + E)e^ t M r pt Gu) = 


D 


2n/i 


2irn 2irn 


rk 


esc 


rk 


rk 2irn\ D 

p + -vr cot ——-1 - E 

2 n rk j r 


2 nfi 

e rk , 


with D and E being the same as (2.30). Then the 1/4 BPS Wilson loop expectation value is 




f 2 \ 


-i/3 Ai' 


Fr 


( 7 r 2 r 2 fc ) 


1/3 ( \ _ _, r 2 —&nr—2 

JV 24 “P 3 r 2 k 


Ai 


( 7 r 2 r 2 fe ) 


-1/3 


l\r k , r 2 —2 \ 
iV “ 24 + ) 


Ai 


+ G r 


(7r 2 r 2 fc ) 


AT - — + 

24 ^ 


3 r 2 k 


Ai 


(7 r 2 r 2 k ) 


~ 1/3 ^-h + w) 


where 


2irn 2irn D 
F r = -— CSC ■ 


2 ^rn 2ttti 

rk rk r ’ r rk CS ° rk 

The 1/2 BPS Wilson loops expectation values are 


(rk 2nn\ D 

-7r cot ——-1- E 

V 2 n rk J r 


= (w!% n ) pt = 


(£),n\ 

1/2 


:n— 1 


27rn 


Ai 


2 r 


esc 


(7 r 2 r 2 k ) 


1/3 / at_ _k _i_ r 2 — 6nr— 2 

JV 24 “C 


3 r 2 k 


rk 


Ai 


( 2 \ —1/3 ( N _ _fc , r 2 -2 

\n 2 r 2 k) l JV 24 3rW 


( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


(6.4) 


(6.5) 
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If the Wilson loops (4.6) and (4.8) are not 1/2 BPS but (4.12) holds, we can still get vacuum expectation 


value of the true 1/2 BPS Wilson loop (4.10) as 


w 


qm,n\ _ 

1/2 /Pt — 


i n -‘c 
2 r 


2 nn 


Ai 


■ CSC 


I 2 \ 1/3 / at k i r 2 —6nr— 2 

\n 2 r 2 k) l iV 24 ' 


3r 2 k 


rk 


Ai 


( 7 T 2 r 2 k ) 


1/3 f at k i r 2 —2\ 
iV “ 24 + 3/W ) 


( 6 . 6 ) 


with c being (5.4). If (4.12) is not true, we only have some matrix model results. 


We expand the above results in the limit N k 1 with r being fixed, and now for the ’t Hooft 
coupling A = N/k there is 

A> 1. (6.7) 


We make expansion of large A and large k. For the 1/4 BPS Wilson loop we have 

127ri + rnr 2 + 24 H n 


/T47( 2 ^+l)’ n \ 

' ”1/4 /Pt 


2nvr V 2 


2 n 2 7r 2 / 1 

3 r 2 k 2 l k A 


2vr 


( 6 . 8 ) 


ti7 r (l2n7ri + 24 + n 2 n 2 — 12r 2 + 36?’n — 48?’ + 24ni7 n ) ^ 


3r 2 k 2 


(*)) 


12a/2A 


+ 0| A 


For the 1/2 BPS Wilson loops we have 

_ /TpM> n \ _ * ^ „mrV2\ 

\ W l/2 /Pt - \ M 1/2 /Pt - - 


2 n 2 7r 2 / 1 

3r 2 fc 2 V k 4 


4n7r 

1 , ti 2 7t 2 + 36rn — 12r 2 + 24 
3 r 2 k 2 


(6.9) 


_l 2 + 




These are in accord with the results in saddle point approach (5.7). Note that for the leading contri 


bution of large k. i.e. the genus zero part, there is no r dependence. 

6.2 1/2 BPS Wilson loops in hook representations 

Now we turn to half BPS Wilson loop^] in hook representations based on |25j, where the density 


matrix for Fermi gas dual to ABJM theory was obtained as (2.35). For the J\f = 4 orbifold ABJM 
theory, we have 

p r = p r . (6.10) 


Similar to the calculation in 1251, we obtain the following result in J\f = 4 orbifold ABJM theory, 


(*)<n 


f(e*) v GC 


\ UU 

/ =Det(l + zp f , r )), 


( 6 . 11 ) 


We have the i/u-loops (4.6 1 and i/^-loops W/T 14.81 with i = 0,1, ■ • • , r— 1. If all of them are half BPS quantum 


rW 


mechanically, we can calculate their vacuum expectation values in the matrix model as shown in this subsection. Due to 
the Z r symmetry of the theory, the results are independent of l. From now on we can omit the index i and subscript 1/2. 
We add subscript r to some quantities of the AT = 4 orbifold ABJM theory to distinguish them from their counterparts 
in ABJM theory. Also the results are the same for the i/ii-loops and 1 / 2 -loops, and so we will not write the same results 


twice. If only a special combination of the 1 / 1 -loops and 1 / 2 -loops (4.101 is half BPS quantum mechanically and (4.121 


holds, the following calculations still apply provided that a constant c (5.41 is added to the result. In the worst condition 


(4.121 does not hold quantum mechanically, and the calculations here are just matrix model results. 
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where 


( 6 . 12 ) 


Pf,r ~ PfP > 


with pf being the same as (2.37). 

The generating function for the half BPS Wilson loop in hook representations (a\b) was given 
in 137] 


1 + (s + t) J2 s a t b W T) 


(a|6) 


(6.13) 


a,6=0 


N 


=Sdet ii±§)=n (1+te " >(1+se " J) 


slJJ (1 — seW)( 1 — te v :i )' 

Therefore, the grand canonical ensemble expectation value of a circular 1/2 BPS Wilson loop W r r a \b) 
in the J\f = 4 orbifold ABJM theory in Euclidean space becomes 


1 + (s + t) ^ s a t b W r ^ a |b)y 


a, 6=0 


\GC det(l + zpfp r x ) 
det(l + zp r ) 


= det 


OO A.rr* _ 1 

a,b Z P 




a,6=0 


oo 1 

A, , Z P 


! + („ + *) 


a, 6=0 


with the states |a) and (6| being defined the same as (2.39). We get the relation 


(W r , {alb) ) GC = (a\f^\b) 

~r— 1 


zp 


i — 1 


+ zp r 

zp' A „ / 1 


= Tr (WIE—rlk>( a l) =Tr (— ft— - -e H \b)(a 

V 1 + Zp r J \ e rH-n _j_ l 

Using Sommerfeld expansion we get 

(W r ^a\b)) GC = 7r3 /i csc(7T d^)m r {p), 

where 

m r (p) = TY(0(// — rJ7)e^|6)(a|) 

= Tr {0(p/r — H)e bI \b){a\) = m(p/r). 

Note that for a circular half BPS Wilson loop in ABJM theory kU( a |6) there is 

(W(a|6)) GC = 7r5 M csc(7T d M )m(/x), 

and then we have 


(6.14) 


(6.15) 


(6.16) 


(6.17) 


(6.18) 

(6.19) 
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Then we can use (2.43) and get 


(W r ,(a\b))% C 


9 Jo(o+l) —16(6+1) jn—l 2vrn 2 ^ 

- r -esc ——e rk 

[a]! [6]! 2 r rk 


( 6 . 20 ) 


where n = a + 6 + 1 is the number of boxes of Young diagram (a\b). The 1/2 BPS Wilson loops 
expectation values in the canonical ensemble are 


{^r,{a\b))pt ~ 


q \a{a+l)-\b(b+l) p-1 27m Ai 




2 r 


CSC 


\ —1/3 


]\t k | r 2 — 6nr —2 
iV “ 24 H 3r 2 fc 


rk 


Ai 


( 7 r 2 r 2 fe ) 


-1/3 


!\T_ k_ , r 2 -2 \ 

- /V 24 + 3/^fc / 


In the large iV limit, the expectation values scale as 


(Wr,(a|&))pt ~ 


^nir\/2\ 


( 6 . 21 ) 


( 6 . 22 ) 


6.3 1/2 BPS Wilson loops in non-hook representations 

Let us consider the half BPS Wilson loops in general representations R = (a\ ■ ■ ■ a s \bi ■ ■ ■ b s ). The 
Giambelli formula states that 


^ r r,(aia 2 ---a s \bib 2 ---b s ) > e 0 det W ri ( ap |i, g ) (c^% e "0- 

PiQ 


As in 1251, we consider the following generating function 


(6.23) 


W r (lV) 


(det(5 P9 + tW r>{ap | 6 ,)(e w , e^)))- 


Similar to computations in 1251, with the dehnitions (2.47) and (2.49) we can get 


(AO = JJ[ d i'i]de t p r (i/i,i/j), 

i 

W A N ) = J^j lld'A'lct (Pf{PQY~ l ){vi,Vj). 


(6.24) 


(6.25) 


Note that the multiplication between boldface variables is understood as matrix multiplication with 


indices p. v and summation being replaced by integration with measures [d/r], [dz/] in eq. (2.47). We 
then have 


E r (z:) = ^2 z N Z r {N) = Det(l + zp r ), 

N =0 
oo 

J2 z N W r (N) = Det(l + zpf(P Q) r ~2 ■ (6.26) 

N=0 


Using the relation 

Det(l + zpfiPQy- 1 ) = Det(l + zp r ) det ( 5 pq + zt(a p |(1 + zp r )~ 1 p r ~ 1 \b q )), 


we can get 


{det{5 pq + tW rXap \ bq) )) 


GC 


det (5 pq + zt(a p \(l + zp r ) 1 p r O 6 ?))- 

p,q 


(6.27) 
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The coefficient of t s in both sides of the above equation gives 


(W r ,(a 1 a 2 -a s |6 1 6 2 -b s) ) GC = det(2:(ap|(l + 2p r ) V 1 |*» 9 »- 


p,q 


(6.28) 


Restricted to hook representation cases, we have 


(Wr,(a| b) ) GG = z{a\(l + zp r ) V 1 \b) 


(6.29) 


So hnally we get 


(^r,(aia 2 ---a s |bib 2 ---& a )} det(W rj ( 0p |j,g)) 


(6.30) 


This shows that the grand canonical ensemble expectation values of the circular 1/2 BPS Wilson loops 
are Giambelli compatible at least in the matrix model sense. 


7 Conclusions and discussions 


In this paper, we have calculated the vacuum expectation values of the circular BPS Wilson loops in 
arbitrary representations in the AT = 4 orbifold ABJM theory. We used both the saddle point approach 


in 1221 and the Fermi gas approach in 23,24 , and the results agree with the available gravity results 
in 


31,32 . It will be quite interesting to study the string/M theory dual of the Wilson loops in the 


higher dimensional representations. 


There are other J\T = 4 SCSM theories 29,40 43 and AT = 3 SCSM theories 44-47 . Recently 


there have been investigations of partition functions of these theories in the Fermi gas approach [48|49|. 
There are also 1/2 BPS Wilson loops in general AT = 4 SCSM theories [321. It would be interesting 
to investigate the vacuum expectation values of supersymmetric Wilson loops in these theories. 

Expectation values of the 1/4 and 1/2 BPS Wilson loops of orbifold ABJM theory in weak coupling 


can be calculated directly in the matrix model, like the ABJM theory case in 17 . Also one can 
calculate the vacuum expectation values of Wilson loops perturbatively using Feynman rules in the 


orbifold ABJM theory in weak coupling, like the ABJM theory case in 12 14,50-52 . It would be 
nice to compare the results of the matrix model to the results of Feynman rules. In fact it was 
proposed in [32] that a perturbative calculation of expectation values of the 1/2 BPS Wilson loops 


using Feynman rules would be helpful in fixing the coefficients in the true 1/2 BPS Wilson loop (4.10). 
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A More general 1/2 BPS Wilson loops in orbifold ABJM theory? 


For two adjacent gauge fields in the quiver diagram and matter fields that couple to them of the AT = 4 


orbifold ABJM theory, one can define two kinds of 1/2 BPS Wilson loops, i.e. the 'i/q-loop 131 32 and 


the ^ 2 -loop 32 . In this appendix we will investigate if there is more general 1/2 BPS Wilson loop 
that preserves the same supersymmetries as the Vh-loop and ^ 2 -loop. 

There is no spacelike BPS Wilson loop in Minkowski spacetime |;39!j. The BPS Wilson loops 
along straight lines in Euclidean space are just the timelike BPS Wilson loops of straight lines along 
Minkowski spacetime after Wick rotation. The circular BPS Wilson loops in Euclidean space can be 
obtained by the conformal transformation of the BPS Wilson loops along infinite straight lines. Also, 
for straight lines the cases of Poincare supersymmetries and conformal supersymmetries are separated 
and very similar. So it is enough to just consider the Poincare supersymmetries of the 1/2 BPS Wilson 
loops along timelike infinite straight lines in Minkowski spacetime. 

We use the conventions in 31,39 . Especially we choose the coordinates x M = (x'/x/x 2 ), and we 
use the gamma matrices 

= (io-W 3 ), (A.l) 

with ( 7 1,2,3 being the Pauli matrices. For the infinite straight line x^ = rd(J, we want to get a 1/2 BPS 
Wilson loop that preserves the Poincare supersymmetries 


70 o u = if? 1 / l 0 e 2i = -hi 2 / 

4~7o = i#K, to = — i 4, 


(A.2) 


A A (2£— i— (2£^) 

with i = 1,2. We only use the gauge fields A/ ; and A/ ; and matter fields that couple to them. 
A general Wilson loop would be of the form 


4/2 = F’exp |^-i J d tZ// 2 (t) ) , 
y(2i+l) j{M) \ 

/f> A< 2 ‘> )' 


At, _ ( -4< 2,+1 > /T> 

-^1/2 “ 


^(2' + d= 4 2 ' +i)iM +x (^rv^i)+ 

A (2C) = A^x" + y (Nifyu- l*l» 


(A.3) 


fi 2t) = V 

f [ 2 2i) = 


—(2£) • . • —(2£) (2£) 
with 77 V ' and r/ 2 ^> being Grassmann even spinors. To make the loop BPS we must find J and g\ 
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that satisfy 153 


M (»i, = i(j f > 9 f> _ ^ 

SAW = i(/f>,?'> - 9 ?' ) /r ) ), (A. 4) 

sfW = p T s! 2() ^ a T si 2<) + iA^w s w _ 

4 /f 1 = Pxsf 1 = a T sf> + iAWgW - i 9 fU< 2 '+‘>. 

Because of the form of f[ 2 ^ and fP\ terms with fields V , ( 2 ^+i) an d should cancel in the 

variation of A^ +l \ Similarly, terms with V^—i) an d ^ should cancel in the variation of A^. 


This forces us to choose 

If we take the ansatz 

we get the ^ 1 -loop with 


Or if we take the ansatz 

we get the t/^-loop with 


M) = N.A = diag(—1,1). 
vP = rt\ 2€) = 77(2£)<Ti, 

M) = NJ = diag(l, 1), 

fjW) a = V(2£)a = (b l)/3, PP=~i- 

4 2i) = V m)S l V(2e) = V(2£)&2i 


M) = N/ = diag(-l,-l), 


(A.5) 

(A.6) 

(A.7) 


(A.8) 


(A.9) 


jf 2£ ) a = P( i, 1), ?7(2£) a = (-i, 1)P, PP = ~ i- 


We wonder if there is a more general 1/2 BPS Wilson loop that preserves the same supersymmetries 
(A.2) as the V’l-loop and ^ 2 -loop, at least classically. One of the consequences of (A.4) is that 


-(21) -(2i)i ,(2t) (21) (2£)7i 

9 1 = 9i ’ 4>\ , 92 = <4 <W 


(A.10) 


with g^ 1 and g^p being Grassmann odd and having no color index or spinor index. From (A.4) and 
the variation of we must have 


700* = -i M'f - iy T,; 2f)9 f') ! , 

5 Vfi 5 , : . J 2 ^ 

t'ul'o — + !y ■ 


(A.ll) 
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Then from (A.2) we have 


e v = _ M ^-\l^r ] } 2£) -g[ 2e)i 


87 r 


- e 2i = -M\e 23 - \l— 

8ir 


k 2 -(2^)i 


r l(2£)9l 


“li ~ AZ jt/ij + y ^Vi <72i > 

a _ A/ rifl , t fk <2i) (2i) 

-t>2i - -M ^2] + Y ^;h2 ff 2 i • 


(A.12) 


Note that 0 1 * and $ 2 * are nonvanishing, general and linearly independent, and similarly 9n and 0 2 * 
are nonvanishing, general and linearly independent. First of all, pj 2£ - ‘ and cannot be vanishing, 
otherwise there would be no solutions for the matrix Aft. Then we must have r 7 ^ = 0 or 7/ 2 2f j = 0, 
as well as r/j 2 ^ = 0 or r/ 2 2 ^ = 0- When 77 ^ = we have ML = — <5j, and then there is 77 ^ = 0. 
This gives the 0 2 -loop. When i] 2 2 i) = 0; we have M l ~ = S'-, and then there is fj^' = 0. This gives the 
■ 01 -loop. 

In summary we have no choices other than the 01 -loop and ■ 02 -loop that satishes the following 
conditions. 


It is constructed by two adjacent gauge fields A 


( 21 + 1 ) 


and in quiver diagrams and fields 


that couple to them in the general form (A.3). 


It preserves the same supersymmetries as the 01 -loop and 0 2 -loop (A.2), at least classically. 


This result can be taken as a small step towards classification of BPS Wilson loops in J\T = 4 SCSM 
theories. 
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